Chiral Phase Transitions around Black Holes 
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In this paper we discuss the possibility that chiral phase transitions, analogous to those of QCD, 
occur in the vicinity of a black hole. If the black hole is surrounded by a gas of strongly interacting 
particles, an inhomogeneous condensate will form. We demonstrate this by explicitly constructing 
self-consistent solutions. 
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According to the theory of quantum fields in curved 
space, black holes radiate energy at a temperature in- 
versely proportional to their mass As the black hole 
evaporates, its temperature rises, and at some point a 
bubble of a high temperature phase surrounding the hori- 
zon may form, if a phase transition occurs. This was 
a particularly interesting phenomena in connection with 
the Higgs model of electroweak symmetry breaking and 
its study requires the inclusion of interactions, being an 
essential feature of the phase transition. A method to 
deal with this situation was proposed, but the indication 
was that, in the Higgs model, the associated high tem- 
perature phase would be too localized around the black 
hole, so that symmetry, effectively, would not be restored 
0. The same problem has been reconsidered by Moss 
taking into account the effect of trapped particles, i.e. 
particles emitted by the black hole and reflected back 
by the walls of the bubble. He indicated that, for some 
class of bag models, the picture may change and lead to a 
transient equilibrium configuration of restored symmetry 
phase, localized around the black hole @ . 

A field in which the similar problem of understand- 
ing the phase structure is nowadays very topical is that 
of QCD at finite temperature and density, in which 
phenomena like chiral symmetry breaking and confine- 
ment/deconfincmcnt transitions are known to take place. 
In this context, the natural way of addressing the prob- 
lem would be to use 1 first principle'' non-perturbative lat- 
tice methods, but already in flat space, and especially 
at high densities, things become prohibitive. In lack of 
a first principle approach, approximating QCD with a 
strongly interacting fermion effective field theories comes 
in handy. The price to pay is that we have to work with a 
non-renormalizable effective theory, but with the bonus 
of dealing with a simpler one that shares many of the 
essential properties of QCD. As a matter of fact, a great 
deal of attention is currently paid on mapping various 
phases on the temperature-density diagram within such 
an effective field theoretical approach, in order to gain 
understanding of the vacuum structure of strongly inter- 
acting matter (See Ref. Q for a recent review). 

The aim of this work is to use the same simplifica- 
tion of degrading QCD to a non-renormalizable, strongly 
interacting fermion effective field theory, and study the 
interplay with black holes. To begin with, we wish to 
consider a little more in detail the issue of phase tran- 
sitions that would break or restore chiral symmetry. In 



the context of strongly interacting fermionic systems, it 
is well known that chiral symmetry breaking takes place, 
and this fact is discussed in terms of the appearance of a 
fermion condensate. One aspect particularly important 
to us is that the ground state is believed to develop inho- 
mogeneous phases when the density becomes large. For 
instance, Refs. @ discussed the issue of chiral symme- 
try breaking and the related condensate formation, and 
mapped the phase diagram for models of the Nambu- 
Jona Lasinio class. The description of Refs. @ indicated 
that the fermion condensate at high densities rcnscmblcs 
a lattice of domain walls. 

In the present case, we are lifting the situation to 
curved space, where new effects kick in. In a constant 
curvature space, the effect of the non-trivial geometry 
is something similar to adding chemical potential. The 
condensate may or may not be spatially homogeneous. 
By contrast, in a black hole spacetime inhomogeneous 
configurations for the condensate are inevitable. To keep 
the situation as simple as possible, let us concentrate on 
the case in which a Schwarzschild black hole of mass m 
is surrounded by strongly interacting fermions in ther- 
mal equilibrium with the asymptotic temperature given 
by Tbh = (8vrm) _1 . Then, the local (Tolman) tempera- 
ture is given by T; oc = Tbh/V7 with / = I — 2m/r. In 
flat space, the strongly interacting fermionic theory has 
a critical temperature, T c , that marks the phase tran- 
sitions of chiral symmetry breaking (in QCD T c ~ 200 
MeV). Therefore it seems evident that in the asymptotic 
region chiral symmetry is restored when Tbh > T c while 
broken for Tbh < T c . When T B h < T c , Ti oc crosses the 
critical temperature at a certain radius. Within this ra- 
dius, the symmetry will be restored. This indicates the 
possibility that a domain wall structure of the condensate 
surrounding the black hole will arise. 

We will now make the above picture quantitative by 
using a strongly interacting fermion effective field theory 
of the Nambu-Jona Lasinio type. The prototype action 
can be written as 



S = I d 4 x^ i #y"V„V> + 



In the above expression ip is a spinor field, A is the cou- 
pling constant, g = |Detg^„| is the determinant of the 
metric tensor and 7 M are the gamma matrices in curved 
space. The number of fermion degrees of freedom (equal 
to the number of flavors x the number of colours) is N 
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and summation over color and flavor indices is under- 
stood. The background spacctime is that of a spherically 
symmetric and asymptotically flat black hole, 



ds 2 



fdt 2 + r X dr 2 + r 2 (d0 2 + 



The formulae we will present below generally apply for 
any function f(r), but the numerical analysis will be car- 
ried out for the Schwarzschild case. 

To analyze the breaking/restoration of chiral symme- 
try, we will use the finite temperature effective action 
in the large- N approximation. The effective action (per 
fermion degree of freedom), T, can be expressed, after 
bosonization, as 

r = - J d^x^-g {^) + Trln (i 7 "V M - a) , 

where the composite operator a = —j^-'ipip was intro- 
duced and the determinant acts both on field and coor- 
dinate spaces. Chiral symmetry is broken dynamically 
when a acquires a non-zero vacuum expectation value 
and then a fermion mass term appears. 

The computation of the effective action can be per- 
formed using the method described in Ref. [H, although 
some modifications are necessary to include the case of 
black holes. Since black hole spacetimes are static but 
not ultrastatic, we rescale the metric (JTJ) so as to be 
ultrastatic, ds 2 = f~ 1 ds 2 . We will use a hat to indi- 
cate the quantities evaluated in this conformally related 
spacetime. After the conformal transformation, one can 
use the method of @ to evaluate the effective action in 
the rescaled spacetime, T, and add a correction term, ST, 
sometimes called cocycle function, to compensate the ef- 
fect of the conformal transformation [3, H| . Assuming the 
condensate to be spherically symmetric, a = cr(r), and 
squaring the Dirac operator, we obtain 
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In the above expression □ is the D'Alembertian 
in the conformally rescaled spacetime and 
a 2 := a 2 + e/ 1/2 er'. The quantity j^W = 
/((n-2)Aln//4-(n-2) 2 (Vln/) 2 /16) is deter- 
mined SO that □ + st^ = /-(«+2)/4 Q j(2-r l )/4 ig 

satisfied, where n is the spacetime dimensions. In 
(J3j) we have used the notation srf = £/^ 4 K Notice that 
stf = R/6. Imposing the periodicity in the Euclidean 
time with the period /3 = 2ir /Tbh, we express T as 



He £ 

e— ± n— — oo 



Trln 



-A 



with A being the Laplacian in the conformally rescaled 
space and ui n := 2ir/f3 (n + 1/2). 



Using zeta rcgularization gives 

f = \ I d 3 x^~g[aO)\ne 2 + C( 



'') ■ (1) where £ is a rcnormalization (length) scale and 
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The quantities £(0) and £'(0) are the analytically contin- 
ued values of £(s) and its derivative to s = 0. The compu- 
tation of the effective action is rather involved, but it can 
be performed in a straightforward manner following the 
method developed in Ref. [Bj]. Here we use a resummed 
form for the heat-trace and retain all terms that contains 
a specified number of spatial differentiations. We invite 
the reader to consult Ref. [f| for details and further ref- 
erences. In the present case, the result is 
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The other term to compute is the cocycle contribution 
that compensates the difference due to the conformal 
transformation to recover the result in the original space- 
time. The cocycle function can be expressed in terms of 
the heat-kernel coefficients associated to the operator & 
in n dimensions: 
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For an operator of the form G = □ + V, the part of the 
heat-kernel coefficient, relevant for our computation, is 



ck%] 



1 



(4tt)^ 



o / ^ X V~9 



v 2 - \rv 



where the dots represent terms that do not depend on V 
or disappear upon integration by parts. In the present 
case, V = a 2 in the original spacetime while V = £/ +f(J 2 
in the conformally rescaled spacetime. Simple computa- 
tions give 
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where A n 



(R^ - fR^)/6 and R.M and 
arc the n-dimcnsional Ricci scalars in the confor- 
mally rescaled and original spacetimes, respectively. Us- 
ing lim„^ 4 dA n /dn = {f' 2 - 2//" + 4//'/r)/24, one ar- 
rives at the expression for the cocycle function. Com- 
bining ([5]) with the above expression gives the effective 
action T for the condensate a. 
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The problem is now reduced to finding extrema of the 
effective action T with respect to the condensate a. Ig- 
noring fourth order derivatives of the condensate allows 
us to express the equation of motion for the condensate 
as a non-linear Schrodinger-like equation of the form 

a" + Skt' + 6 2 (t' 2 + = , (5) 

where the coefficients Si and are functions of a but 
independent of its derivatives. The explicit expressions 
are rather long and will not be reported here. 

Before finding the explicit solution for the condensate, 
we will discuss the critical temperature in the asymp- 
totic region r — > oo. Denoting the minus of the action 
with a 1 = as the potential U(a), the derivative of the 
asymptotic value can be computed exactly as 

3cr (4Acr(4w_i (er) + /foln (o/i) - 2\(3a 2 + /3) 
" as = 2A/?(-4/?(7OTi( ( 7) - 6w (a) + 3 In (a/£) - 2) 

The critical temperature is determined by the equation 
d 2 U as (a) = 0. Thus, expanding the Bessel functions con- 
tained in sj„ for small a, performing exactly the sum over 
n, and finally solving a trivial algebraic equation, one ar- 
rives at T c = yJHX^ 1 / 2 . The thermodynamic potential 
obtained by numerically integrating d a U with respect to 
<7 is shown in Fig. [T] 




Figure 1. The figure illustrates how, asymptotically, the 
potential U a3 (o') changes and symmetry gets restored as 
temperature increases (The top (red) curve corresponds to 
Tbh /T c = 1.75, while the bottom (orange) curve corresponds 
to Tbh /T c = 0.03. Fhe second curve from top (blue) corre- 
sponds to T/T c = 1.). We set £ = 10 6 and A = 10~ 2 . 

Computing the thermodynamic potential locally will 
provide further insight on the form of the condensate. In 
fact, such a computation shows that starting from a set 
of parameters for which asymptotically the potential has 
a non vanishing minima, as we move towards the hole, 
the minima of the potential will gradually shift towards 
a configuration with vanishing a . We confirm the above 
picture by solving Eq. §5§ for the condensate with regu- 
lar boundary conditions at the horizon. Solving Eq. ([5]) 
can be handled by standard numerical techniques, but it 
requires some caution. First of all, we notice that the co- 
efficients of Eq. ^ for a depend on infinite summations 



over Bessel functions, whose argument is proportional to 
the condensate. When the value of the condensate is 
not small, these sums can be truncated due to the ex- 
ponential fall-off of the Bessel functions. However, when 
the condensate is small, fully resummed expressions have 
to be used. Once we expand the Bessel functions for 
small values of their arguments, we can perform the full 
resummation over n. In the region r < where a is 
small up to a value r* , we integrate Eq. ([5]) using this re- 
summed form, and then we switch to the truncated form 
for r > r*, matching the value of a and its derivative at 
the junction. The boundary conditions in the vicinity of 
the horizon and in the asymptotically far region are set 
by requiring that the condensate is at a minima of the 
potential. 
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Figure 2. The figure illustrates the condensate profile 
found by solving Eq. ©, for four indicative values of 
the black hole temperature (Left to right: Tbh /T c = 
0.50 (Blue), 0.54 (Green), 0.58 (Red), 0.61 (Black)). The val- 
ues of the other parameters are set to I — 10 3 , A = 10 -2 . As 
we increase the black hole temperature, the region of restored 
symmetry phase expands, and the bubble becomes larger and 
thicker. The asymptotic value of the condensate becomes 
smaller as the asymptotic temperatures increases, and tends 
to zero for T — >• T c . The small box superposed illustrates 
for the rightmost curve (T/T c = 0.61), the corrected solution 
(red, dashed) when fourth order derivative terms are included. 
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Figure 3. The figure illustrates the condensate profile for 
values of the black hole temperature much smaller than 
T c . The curves refer to (left to right): T B h/T c = 
0.34 (Red), 0.35 (Orange), 0.37 (Yellow). 
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We present the results for the condensate profile in 
Fig. [3] for sample values of the parameters. The kink- 
type configurations of Fig. [3] are bubbles that separate a 
region of restored symmetry near the black hole from a 
region of broken symmetry surrounding it. The size of 
the bubble can also be easily estimated by equating the 
local temperature to the critical temperature as 

nubble ~ r B / (l - T% H /T%) , 

which approximately agrees with the numerical results. 

Higher order corrections can be treated systematically 
in the present scheme by iteratively including higher or- 
der terms in the heat-kernel expansion. We expand a 
around the solution already obtained in the lower order 
approximation, a, as a — a + da. Substituting this form 
into the equation of motion, and suppressing all deriva- 
tives higher than three acting on 5a, we obtain a second 
order differential equation for 8a with a source term. We 
carried out this iteration to the forth order and verified 
that higher order corrections only produce small distor- 
tions compared with the solution truncated at the second 
order. In fact, higher order terms become less and less 
relevant as the black hole temperature gets closer to the 
critical one, due to the fact that the kink becomes in- 
creasingly thicker. In Fig.[3]we have shown the perturbed 
solution superposed to the one truncated at second order 
for an illustrative set of the parameters. 

The computation presented here considers a situation 
of thermal equilibrium, and thus we used the Hartle- 
Hawking vacuum state. To describe evaporating black 
holes, we have to use the Unruh vacuum state. The 
outgoing flux in the Unruh vacuum state is diluted at 
infinity, effectively lowering the asymptotic temperature. 
This suggests that the formation of a bubble can occur 
also for black holes with temperature higher than T c . 
This case is technically more challenging and work in 
this direction is in progress. 

The present discussion may be of relevance in the con- 
text of primordial black holes. Recent attention has been 
drawn to the possibility of a chromosphere formation 
around the black hole, and Ref. Q, contrary to the orig- 
inal claim of Ref. [l(| , suggested that this does not hap- 
pen. We should stress that the phase transition across 
the kink discussed here is a phenomena largely different 
from the chromosphere formation. For the appearance of 



the symmetry restored phase around a black hole, scat- 
tering or reflection on the phase boundary of particles are 
not so essential. In the context discussed in this paper, 
the gradient of the effective local temperature caused by 
redshift plays an essential role. In the case of the Un- 
ruh vacuum, as mentioned above, geometrical dilution 
of particles assists this tendency. It would be definitely 
interesting to discuss whether the appearance of the sym- 
metry restored phase around a black hole may change the 
discussion of Ref. H, [n| ■ 

Several other interesting generalizations of the present 
work include the case of charged/rotating black holes. 
Analyzing the case of higher dimensional black holes 
would also be interesting in view of the possibility that 
microscopic black holes may form at the LHC. 

Here we considered the simplest class of models and 
completely ignored the role of gauge degrees of free- 
dom or other types of condensates, such as pseudo-scalar 
ones. Improving the description in this direction is es- 
sential and may help us to gain insight into possible con- 
finement /deconfincmcnt transitions. Again, the analogy 
with QCD turns out to be of use and, in the context of 
the effective theory approach adopted here, one natural 
possibility is to couple the model to the Polyakov loop 
fill ] - While adding a pseudo-scalar condensate does not 
produce significant changes, the inclusion of gauge fields 
requires some additional efforts. Work in this direction 
is in progress. 

From this point on, speculation would bring us too 
far, but certainly there are various interesting problems 
to consider in relation to the discussion presented here. 
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